arXiv: 1505.06823vl [cond-mat.mes-hall] 26 May 2015 


High-Tc dP ferromagnetism in a doped Mott insulator: the case of 
hydrogenated epitaxial graphene on SiC(OOOl) 

Pengcheng Chen 

Department of Physics and State Key Laboratory of Low-Dimensional Quantum Physics, 
Tsinghua University, Beijing lOOOSf, People’s Republic of China and 
Collaborative Innovation Center of Quantum Matter, 

Tsinghua University, Beijing 100084, People’s Republic of China 

Yuanchang Li* 

National Center for Nanoscience and Technology, 

Beijing 100190, Peoples Republic of China 

Wenhui Duan^ 

Department of Physics and State Key Laboratory of Low-Dimensional Quantum Physics, 
Tsinghua University, Beijing 100084, People’s Republic of China 
Collaborative Innovation Center of Quantum Matter, 

Tsinghua University, Beijing 100084, People’s Republic of China and 
Institute for Advanced Study, Tsinghua University, 

Beijing 100084, People’s Republic of China 
(Dated: May 27, 2015) 


1 


Abstract 

We show that the dP ferromagnetism with high Curie temperature (Tc) can be achieved in the 
electron doped hydrogenated epitaxial graphene on some certain SiC substrates through first- 
principles calculations. The pristine systems are found to be a Mott insulator independent of 
SiC polytypes {2H, AH or 6i7) which, however, play a significant role in the modulation of mag¬ 
netic interaction. Carrier doping enhances the ferromagnetic coupling due to the double exchange 
mechanism and thus realizes the phase transition from antiferromagnetism to ferromagnetism. A 
Tc of around 400 K is predicted on the 2i7-SiC. We employ a non-degenerate Hubbard model to 
demonstrate how the SiC affects the interfacial magnetism in intra-atomic Coulomb repulsion and 
intersite hopping interactions. 
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ferromagnet particularly with high Curie temperature (Tc) is a long-sought goal for 
its nontrivial magnetic mechanism and long spin relaxation time which would promise a 
better performance in future spintronic devices[l-3]. In low-dimensional systems, the elec¬ 
tron kinetic energy is decreased with the reduction of atom coordination number at the 
surface or interface having localized defect states, resulting in an enhanced Coulomb in¬ 
teractions/bandwidth ratio. This usually gives rise to a large density of states around the 
Fermi energy and thus the emergence of spin-splitting according to the Stoner criteria[4, 5]. 
Following such a paradigm, ferromagnets have been predicted based on 2p-electrons of 
carbon, oxygen and nitrogen at the interface, surface, edge or defects[6-13]. Basically, cur¬ 
rent semiconducting industry is established on the Si element. Thus the integration of d^ 
ferromagnets with Si-based systems is also highly desired. 

Recently, room temperature d^ ferromagnetism was reported in hydrogenated epitaxial 
graphene on SiC(0001)[14, 15], yet the physical origin is still ambiguous, especially about 
the interpretations of local moments. Some researchers attributed it to the formation of 
unpaired electrons in graphene[14] while others ascribed it to the interfacial Si dangling 
bond (DB) states[15]. This issue is also crucial from a practical point of view: if the former 
is true, an additional protection process has to be implemented to retain the unpaired 
electrons in the presence of air; otherwise, the graphene layer itself would offer an inherent 
protection. On the other hand, it is a bit surprising that little attention is paid to the 
SiC substrate used although its polarization effect was reported to be stacking sensitive, 
which naturally affects the local electrostatic environments at the interface, e.g., the on¬ 
site Coulomb repulsion[16, 17]. In fact, different SiC polytypes (namely, 4Fr-SiC[14] and 
6i7-SiC[15]) were employed in the experiments for the epitaxial growth of graphene. 

In this work, we report the engineering of d^ ferromagnetism in the hydrogenated epi¬ 
taxial graphene on different polytypes of SiC through hrst-principles calculations, where the 
moment is associated with the interfacial Si-DB state. Importantly, we hnd a new route to 
tune the magnetic interactions between the local moments, namely, the stacking sequence 
of SiC bilayer. The pristine system, regardless of SiC stacking, is a Mott insulator and 
possesses an antiferromagnetic ground state but the magnetic coupling strength depends 
upon the stacking heavily. When the carrier density increases, the ferromagnetic coupling 
is enhanced because of the double exchange mechanism and an estimated about 400 K is 
achieved at an electron density of 4xl0^^/cm^ on 2i7-SiC. We employ the non-degenerate 
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Hubbard model to demonstrate how the SiC stacking affects the magnetic properties of the 
system. Our work not only unravels the origin of ferromagnetism in hydrogenated epitaxial 
graphene but also reveals a feasible system that is integratable with the current semicon¬ 
ducting industry for the applications in spintronics. 

First-principles calculations are performed using the projector augmented wave[18] meth¬ 
ods within the local spin density approximation (LSDA)[19] as implemented in the Vienna 
ab initio simulation package (VASP)[20]. The slab model is used with a y/S x \/3R30° 
supercell for all 2H-, AH- and GhT-SiC (0001) surfaces, which can accommodate a 2x2 
graphene cell. The SiC substrate is modeled by justihed four (2if, AH) or six (6iJ) SiC 
bilayers with the H-terminated bottom surface. The SiC bilayers in the lower half are hxed 
at their respective bulk positions. All the other atomic positions are fully optimized using 
the conjugate gradient algorithm without any symmetry constraint until the residual forces 
become smaller than 0.01 eV/A. A Monkhorst-Pack /c-point mesh of 12x12x1 was used for 
all Brillouin-zone integrations. 

Local moments.-Yoi epitaxial graphene grown on SiC(OOOl), the first layer is so called 
“buffer layer”, where a part of carbon atoms chemically bond to the Si atoms underneath. 
This results in the atomic corrugation[21-24], which is typically illustrated in Fig. 1(a) by 
taking 6iP-SiC as an example. It can be seen that two thirds of surface Si atoms bond to 
the carbon in the buffer layer while the remaining one third have the unsaturated DBs. The 
bonding between carbon and surface Si atoms violates the linearly dispersive band structure 
of graphene, leading to the metallic states contributed by Si^ and C^ around the Fermi level, 
as shown in Fig. 1(b). It is the hybridization between Si-DB and C^ states that maintains 
the spin-degenerate states at the interface. 

We then employ the hydrogenation approach to passivate the unsaturated C„ states 
and decouple Si-DB and C„ states without breaking the interface structure[25, 26]. Our 
energetic calculations show that the hydrogenation of yields an energy gain of 0.76 eV/H 
as compared to an energy cost of 0.95 eV/H for Cf,. Therefore, we only consider the situation 
with the six fully hydrogenated as shown in Fig. 1(c) hereafter. Figure 1(d) plots the 
corresponding band structure without taking spin-polarization into account. It can be seen 
that the C„ states are decoupled from Si„ states and shift to lower energy. While around 
the Fermi level, there appears a nearly flat band dominantly contributed by Si^ DB state. 
According to the Stoner criterion[4, 5], spontaneous spin-polarization might occur given the 
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FIG. 1. (Color online) (a) Side and top views of the first graphene layer on 6i/-SiC(0001), where 
the snbscripts of “b” and “n” respectively refer to the tetra- and triple-coordinated atoms, and 
(b) the corresponding band structure, (c) Geometries after hydrogenation and the corresponding 
band structures without (d) and with (e) spin-polarization. Black rhombuses denote the surface 
\/3 X \/3 R30° supercell in (a) and (c). The bands highlighted in red and green represents the 
states mainly contributed from Si^ (Si-DB)and atoms, respectively. In (e), black solid and blue 
dashed lines correspond to the majority and minority spins, respectively. The Fermi level is set to 
zero. 

large density of states at the Fermi level. Indeed, the spin-polarized state is energetically 
more favorable by 75 meV than the non-spin-polarized one. Fignre 1(e) shows the spin- 
resolved band strnctnre, in which the originally spin-degenerate Si„ DB bands fnlly split 
into an occnpied (majority) and an nnoccnpied (minority) band with an energy gap of 
~0.17 eV, hence contribnting to 1 hb local magnetic moment. 
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It is worth to note that although the above results are obtained for 6iJ-SiC, the same trend 
holds for either 2H- or 4if-SiC according to our calculations, and is also independent of the 
SiC bilayer stacking sequence. An important difference is the magnitude of the bandwidth 
and spin-spliting of Si-DB band, which signihcantly affects the magnetic coupling between 
local moments. Below, we will discuss this in detail. 

Long-range magnetic order.-For practical applications, the existence of local moment is 
not enough and the long-range ferromagnetic order must be present. In fact, the emergence 
of charge carriers is unavoidable in our system due to the doping effect from SiC substrate[16, 
17], which will play a signihcant role in mediating the magnetic interactions as intensively 
studied in dilute magnetic semiconductors[27] and ferromagnets[8, 28]. Moreover, it is 
possible to precisely control the carrier type and concentration by applying gate voltage. 
Hence we next investigate the possible long-range magnetic order in our systems while 
including the effect of charge doping. 

To investigate the magnetic coupling of the local magnetic moments, we have explored 
two collinear magnetic structures. One is ferromagnetic (FM), i.e., all the local moments 
have a parallel alignment as shown in Fig. 2(a). The other is collinearly antiferromagnetic 
(ColAFM) with the row-wise antiparallel alignment as shown in Fig. 2(b). The magnetic 
interaction thus can be conveniently mimicked by mapping their respective total energies 
(FM and ColAFM) to the Heisenberg model. Note that the ColAFM conhguration should 
not be realistic on a triangular lattice due to the frustration[29, 30]. However, our test 
calculation on a more realistic coplanar Neel (120°) conhguration (See Fig. 2(c)) gives a 
negligible energy lowering, less than 0.3 meV/Si-DB. Given the huge computational cost 
for noncollinear calculations, only the ColAFM situation is considered to simulate the AFM 
order. 

Figure 3(a) summarizes the energy difference (Ai7) between the AFM and FM states as 
a function of carrier density (p) in hydrogenated epitaxial graphene on different hexagonal 
polytypes of SiC (277, 477 and 677). For 477-SiC, two kinds of stacking sequence are calcu¬ 
lated as schematically shown in Fig. 3(b). Note that the charge doping range is comparable 
with the intrinsic doping level of epitaxial graphene [31] and also experimentally accessible in 
hydrogenated epitaxial graphene [15]. The largest doping situation corresponds to a carrier 
density of 0.1 e/Si-DB (~4x 10^^/cm^). Several features of (A77) curves are observed in Fig. 
3(a). 
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FIG. 2. (Color online) Schematics of the calculated magnetic configurations; (a) Ferromagnetic 
(FM), (b) Collinear Antiferromagnetic (ColAFM) with the row-wise antiparallel structure, and (c) 
Coplanar noncollinear Neel (120°) structure. Dash rhombuses are the unit cells. 

Firstly, the AFM state has the lower energy under the ideal situation without any carrier 
doping independent of SiC polytypes. This can be understood in terms of Fig. 3(c). The 
Si-DB state is fully occupied in the majority spin and the virtual hopping is allowed when 
the nearest-neighboring local moments are in antiparallel alignment but forbidden for the 
parallel alignment conhguration. Therefore, the ground state is AFM at p = 0. 

Secondly, as |p| increases, the system becomes “more” ferromagnetic either for n-type 
(p>0) or p-type doping (p<0). The corresponding physics is illustrated in Figs. 3(d) and 
3(e), revealing a typical double-exchange mechanism. Only when all the spins orient in 
parallel, the doped electrons/holes can move freely from one site to another, consequently 
reducing the system energy. Note that it is the competition between the double exchange 
and antiferromagnetic virtual hopping that determines the system ground state. Thus there 
exists a critical p, over which the AE becomes positive and the system undergoes a phase 
transition from the AFM to FM state. The asymmetric behavior between electron and hole 
indicates the different intra-site hopping interaction of varied carrier types. 

Finally, there exists an interesting linear dependence of AE on p regardless of SiC poly¬ 
type, stacking sequence and carrier type. But different SiC substrates exhibit the different 
slops and intercept (p = 0), showing their distinct effects on magnetic interaction. For ex¬ 
ample, the parallel coupling is enhanced more rapidly in 2H- and 4if2-SiC as the p increases 
in comparison with AHi- and bFT-SiC. At p = 0, the AE increases in the order of 4Fr2-, QH-, 
2H and 4ifi-SiC. These clearly show that SiC substrates, not only the polytypes but also 
the stacking sequence, can be used as an efficient way to modulate the long-range magnetic 
interaction in such a system. 
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FIG. 3. (Color online) (a) Energy difference (AE) between the AFM and FM states as a function 
of carrier density (p) in hydrogenated graphene on different SiC substrates, (b) Schematics of the 
SiC bilayer stacking sequence. The shaded area denotes graphene layer, (c)-(e) Schematics of the 
exchange mechanism without and with carrier doping. For the neutral system, virtual hopping 
yields an antiparallel alignment of nearest-neighbor local moments (c). In contrast, hopping of 
doped electrons (d) or holes (e) mediates a parallel alignment of local moments. 


In fact, the system can be described by the classical non-degenerate Hubbard model 


H = -k UEni^riii, 


( 1 ) 


where the hrst term denotes the hopping (tij) interaction of local electrons from site j to i, 
and the second term denotes the on-site Coulomb repulsion (U) with rii^ = Without 

carrier doping, the AFM exchange arises from virtual hopping of an electron with spin a to 
a neighboring site occupied with an electron of spin -a(see Fig. 3(c)). The exchange energy 
is —2t^/U with the hopping matrix element tij = t. While in the system with a certain 
amount of carriers, their hopping from site to site with the matrix element fgff would lead 
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to the formation of the energy band: 

e{k) = -2fefr[cos(y + + cos(y - + cos(A;a;)] (2) 

for the FM triangular lattice (Fig. 2(a)) and 

e{k) = -2fefr cos{k^) (3) 

for the ColAFM configuration (Fig. 2(b)). Note that the carrier hopping is only allowed 
between the sites with parallel spins. When the |p| is relatively small, just the lowest energy 
level would be occupied, thus approximately yielding the energy gains of = —6fefr|p| 
and i^coiAF = —2fefr|p|, respectively. Further taking the virtual hopping contribution into 
account, we then have 

^E = -2e/U + At,^\p\. (4) 

It clearly reveals the linear dependence of Ai? on p in Eq. (4). At p = 0, the AE 
corresponds to -2t^/U. In fact, t and U can be estimated as the bandwidth and spin 
splitting of Si-DB band as shown in Figs. 1(d) and 1(e), respectively. The derived t {U) 
values are 26.6, 37.0, 47.6, and 63.0 meV (688, 670, 520, and 562 meV) for 4i7i-, 2H-, 
6H- and 4i72-SiC. It yields the value of 2.76, 5.52, 10.1, and 17.2 meV for 2t^/U, in good 
agreement with the data obtained from DFT calculations as shown in Fig. 3(a). Note that 
the ratio of is at least 9 to 1, and therefore, the pristine hydrogenated system should 
be characterized as a Mott insulator. 

Doping charge carrier in Mott insulator will signihcantly alter its magnetic order and 
transport properties, leading to unconventional quantum phenomena such as magnetic phase 
transition[32] and high-temperature superconductivity[33]. The slope fes in Eq- (4) repre¬ 
sents the energy gain due to hopping interaction of carriers. By fitting the data in Fig. 3(a) 
linearly, we obtain three typical values of tes, namely, ~2 meV, ~25 meV and ~60 meV. The 
value for hole doping is nearly universal (~2 meV) while there exist two typical values for 
electron doping on 4i7i- and 6H- (~25 meV) vs. 4772- and 277-SiC (~60 meV). The distinct 
behaviors of electron and hole doping are generally found in doped Mott insulators, however, 
its physical origin is still under debate. [33] In this case, the ferromagnetic exchange energy 
between two Si-DB spins can be evaluated by the form of double exchange interaction[34]: 

where thyb is a hybridization matrix element between the local Si-DB states and 
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surrounding charge densities which connect the inter-site Si-DB states, and A is the en¬ 
ergy difference between Fermi level and energy level of surrounding charge densities. The 
electron-hole asymmetry might he on the different Fermi level and the hybridization inter¬ 
actions. Anyway, the exact dependence of related parameters on the SiC stacking deserves 
to be further explored. 

It is found that except the SiC polytype, the Ai? also depends on the SiC bilayer stacking 
sequence as demonstrated in the comparison between AHi and AH 2 cases (See Fig. 3(a)). The 
trend is also observed on 6hf-SiC (not shown here). According to the previous studies[16, 17], 
the different stacking faults (i.e., an altered stacking sequence with respect to that in ideal 
cubic SiC) cause an entirely different local electrostatic environments at the SiC surface. 
This might be the underlying physics responsible for the modulation of interfacial magnetic 
interaction via SiC substrate. 

Discussion.-\i should be emphasized that the recently reported room-temperature dP 
ferromagnetism in hydrogenated epitaxial graphene[15] is essentially different. Therein, an 
additional graphene overlayer, locating on the first graphene layer, is necessary for the long- 
range FM order, which, however, does not exist in our system. Nevertheless, our work can 
still provide some useful insights for the understanding of the experimental observations, 
e.g., the origin of local moments. More specially, Giesbers et a/. [15] unambiguously clarihed 
the absence of FM signal in the samples with only the hrst graphene layer hydrogenated 
on 6i7-SiC. Indeed, our results shown in Fig. 3(a) indicate that the AFM state is more 
favorable in energy than the FM state at the experimental doping level (0.9 /is/hexagon 
projected area). This is well consistent with the experimental measurements[15]. Whereas, 
we predict that the FM order can be achieved by the change of SiC substrate. For exam¬ 
ple, the Ai7 increases up to 22.6 meV at the same doping level on 2i7-SiC. To estimate 

Tc, we employ the effective exchange Hamiltonian, H = — Jefr X] ^iSj, where Jeff is the 

(ti) 

effective nearest-neighbor exchange coupling and can be extracted from Ai7, the (qj) de¬ 
notes nearest-neighbor sites and Si is the spin value at the site i. Under the mean field 
approximation, we can calculate the Tc through the formula Tc = -^zS{S + 1) Jeff, where 
z = 6 for the triangular lattice. S' = |, and ks is the Boltzmann constant. It gives a Tc 
of 394 K. Combining the high Tc and inherent protection by atop graphene layer, the 
ferromagnetism on SiC should be a promising system for future spintronic applications. 

In conclusion, we have proposed a new route to engineer ferromagnetism with high Tc 


10 



via the modulation of SiC bilayer stacking sequence using the first-principles calculations. 
It is found that hydrogenation always recovers the spin-splitting of Si-DB state while its 
local environment strongly depends on the SiC bilayer stacking sequence. By electron dop¬ 
ing, the long-range ferromagnetism can be readily established. These new results call for 
experimental verihcation. 
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Nos. 2011CB921901 and 2011CB606405), and the National Natural Science Foundation of 
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